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ACOUSTIC  WAVE  SCATTERING  FROM  A PROGRESSIVE 
INTERNAL  WAVE  FIELD 


P.  A.  Barakos 


ABSTRACT 

'^The  scattering  of  sound  from  an  internal  wave  field  is  considered. 

The  acoustic  index  of  refraction  depends  on  the  density  of  sea  water. 
Hence,  spatial  and  temporal  density  variations  arising  from  the  presence 
of  internal  waves  will  produce  a perturbation  in  the  index  of  refraction. 
Therefore,  an  acoustic  wave  transmitted  through  a stratitied  fluid  in 
which  internal  waves  are  present,  will  be  modulated  by  the  internal  wave 
field  and  scattering  and  refraction  will  occur.  It  is  shown  that  the  fre- 
quency spectrum  of  the  scattered  acoustic  field  is  comprised  of  the 
fundamental  acoustic  frequency  and  two  side  frequencies.  Equations 
describing  the  phenomenon  are  given. 
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INTRODUCTION 
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Propagation  of  sound  in  surface  and  subsurface  ducts  has  received 
extensive  treatment  in  the  literature  (1,  2,  3,  4,  5,  6).  Recently, 
Schulkin  (7)  has  summarized  our  knowledge  of  this  important  branch 
of  Ocean  Acoustics.  Perfect  sound  channels  are  strictly  horizontal 
and  their  effects  upon  acoustic  wave  propagation  simply  reduce  to 
certain  anomalies  in  wave  refraction.  In  the  real  ocean,  however, 
imperfect  sound  channels  with  space  - and  time  - varying  boundaries 
are  nearly  always  a universal  phenomenon.  The  imperfections  are 
caused  by  sea  roughness,  bubble  layers,  variations  in  density  grad- 
ients, absorption,  turbulence,  water  currents , current  gradients , 
temperature  microstructure,  biological  effects  and  other  disturbances 
in  the  thermocline  mostly  due  to  internal  waves  and  turbulence.  To 
elucidate  and  overcome  these  imperfections  in  sound  propagation  is 
one  of  the  principal  problems  of  ocean  acoustics. 

Moreover,  the  literature  on  the  theory  of  scattering  of  acoustic 
radiation  from  rough  surfaces  is  very  extensive  (8,  9)  and  the  problem 
of  acoustic  wave  propagation  in  imperfect  channels  has  not  been  com- 
pletely solved.  Clay  (10),  Bucker  (11)  and  Bucker  and  Morris  (12) 
have  considered  some  of  the  effects  of  surface  roughness  on  the  normal- 
mode description  of  duct  propagation  and  find,  in  general,  that  sea 
surface  roughness  contributes  to  the  attenuation  of  the  normal  modes. 

The  related  problem  of  amplitude  and  phase  modulation  of  an 
acoustic  wave  field  from  a periodic  free  surface  has  been  investigated 
experimentally  by  Roderick  (13)  and  Roderick  and  Cron  (14)  in  a model 
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tank.  They  show  that  the  frequency  spectrum  of  the  forward  scattered 
- modulated  signal  consists  >i  a discrete  frequency  spectrum  that  is  cen- 

tered about  the  carrier,  and  a continuous  spectrum  that  is  positioned 
symmetrically  about  the  carrier. 

Similarly,  for  back  scattering  from  a free  surface  in  motion.  Marsh 

. (15)  employing  linear  theory  shows  that  the  reverberation  spectra  are 

narrow  and  centered  at  frequencies  u)0  ± fiwith  D = (2  u>0  g cos  0o/c)^  *n 
which  g is  the  acceleration  due  to  gravity,  0O  is  the  backscattering  angle 

of  the  reradiation,  u>0  is  the  angular  frequency  of  the  carrier  and  c is 

its  phase  velocity.  Lieberman  (16)  confirms  this  prediction  for  sound 


in  air  at  a carrier  frequency  of  78  kHz  over  a randomly  agitated  water 
surface.  For  capillary  waves  having  wavelengths  of  0.  19  to  0.  36  cm 
(for  capillary  waves  v = (2ir  y / A p)^,  where  v is  the  capillary  phase 
velocity,  and  A is  their  wavelength,  p is  the  density  and  y is  the  surface 
tension  of  the  water)  he  predicts  a Doppler  shift  of  167  Hz  with  a dis- 
persion of  30  Hz  about  the  shift.  While  Mellen  (17)  employing  carrier 
frequencies  of  85  kHz  and  1 . 4 MHz  underwater  at  the  Thames  River 
estuary  with  wind  roughened  surfaces  failed  to  observe  the  predicted 
Doppler  spectrum  of  surface  backs catter.  He  hypothesizes  that  the 
apparent  failure  of  first  order  linear  theory  to  predict  the  Doppler 
spectrum  of  surface  backscatter  suggests  that  nonlinear  surface  effects 
may  be  dominant  at  these  wavelengths.  However,  one  should  note  that 
the  first  order  linear  theory  successfully  predicts  the  Doppler  spectrum 
for  both  back-  and  forward  scattering  for  results  obtained  in  model  tank 
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experiments  and  more  recently  for  results  obtained  in  the  ocean,  Roderick 
and  Cron  (14),  at  750  Kz  and  1. 5 kHz.  Nevertheless,  it  is  possible  that 
nonlinear  effects  may  be  important. 

Herein,  consideration  will  be  principally  confined  to  the  lower  boundary 
of  the  imperfect  duct  along  which  internal  gravity  waves  may  be  propagated. 
The  internal  wave  field  will  modulate  the  amplitude  of  the  acoustic  wave 
field  traversing  the  perturbed  interface  so  that  the  spectrum  of  the  scattered 
acoustic  field  contains,  along  with  the  fundamental  acoustic  frequency,  two 
additional  side  frequencies. 

The  significance  of  the  present  investigation  is  that  the  amplitude  and 
phase  modulation  of  the  acoustic  wave  field  by  internal  waves  will  result 
in  sidebands  similar  to  those  observed  by  Roderick  for  forward  scattering 
and  Lieberman  for  backscattering  but  lower  in  frequency.  Consequently 
in  a medium  in  which  both  free  surface  gravity  waves  and  internal  waves 
are  present,  the  resulting  spectrum  of  the  scattered  frequency  will  be 
comprised  of  the  fundamental  acoustic  frequency  and  four  additional  side 
frequencies;  the  outer  two  resulting  from  surface  gravity  waves  and  the 
inner  two  from  internal  waves.  In  the  results  of  project  MIMI,  Steinberg 
(18),  there  seems  to  be  some  experimental  evidence,  Williams  and  De- 
fer rari  (19), to  support  this  hypothesis. 

To  obtain  an  order  of  magnitude  of  the  Doppler  shift  for  internal 
waves,  we  consider  the  two  layer  model.  For  this  model,  internal  waves 
have  their  maximum  amplitude  at  the  interface.  If  the  wavelength  (A)  is 
small  with  respect  to  the  thickness  (H)  of  the  upper  layer,  the  phase 
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velocity  of  the  wave  is  given  by  c = (gAAp/2irp)E  where  g is  the  accelera- 
tion due  to  gravity,  Ap  is  the  density  difference  between  the  juxtaposed 
layers  and  p is  the  mean  density.  On  the  other  hand  if  the  thickness  of 
the  upper  layer  is  small  with  respect  to  the  wavelength,  the  phase  velocity 
is  given  by  c = (gHAp/p)^.  In  the  open  ocean  Ap  = 0.  002  gm/cm3  is  an 
upper  limit  while  for  other  special  situations  values  as  high  as  0.  1 
gm/cm3.  Degens  and  Ross  (20),  may  be  encountered.  The  distribution 
of  wavelengths  (A)  for  internal  waves  varies  widely  from  a few  tens  of 
meters  in  coastal  waters  and  semi-enclosed  basins  to  a few  hundreds 
of  meters  to  many  kilometers  for  the  continental  shelves  and  the  open 
oceans.  From  these  equations,  however,  one  obtains  an  approximate 
value  of  150  cm/sec  for  the  phase  velocity  for  internal  waves.  Hence, 
for  acoustic  radiation  with  carrier  frequencies  in  the  audible  part  of  the 
spectrum,  Doppler  shifts  of  a few  tens  of  Hertz,  at  most,  are  expected. 
MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 

The  procedure  followed  herein  is  exactly  identical  to  that  employed 
in  the  diffraction  of  electromagnetic  waves  by  ultrasonic  waves,  Morse 
and  Ingard  (21)  and  Born  and  Wolf  (22). 

To  investigate  the  scattering  of  an  acoustic  wave  by  an  internal  wave 
field,  we  start  with  the  wave  equation  for  the  acoustic  field 


where  A is  the  two-dimensional  Laplacion,  \|i  is  the  space-  and  time- 
dependent  acoustic  velocity  potential  and  c is  the  speed  of  sound  in  the 
ocean. 
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We  write  for  the  acoustic  index  of  refraction  n = c0/c,  where  c is  the 
speed  of  sound  at  some  reference  depth  of  the  water  and  c is  the  speed  of 
sound  at  some  other  depth.  We  express  the  index  of  refraction  as  the  sum 
of  an  unperturbed  value  n0  and  a perturbed  value  nx  caused  by  the  internal 
waves.  Since  nx/n0  « 1,  then  n2  = n§  + 2n0nx  and  the  wave  equation  is 
expressed  approximately  as 


A2  if  - 2n*/n°  2ll  - 0. 
c2  At2 


If  we  introduce  the  relative  variation  in  the  temperature  which  is 
related  to  the  index  of  refraction  approximately  as  n ~ 1/Tso  that 
26n/n  = 2T/6T  then  Eqn.  (2)  becomes 


a 1 d2i|r  (2T/61)  d2i|i 

A ¥ “ — 


c2  St2  c2  St2 


We  assume  the  coordinate  geometry  of  the  internal  and  acoustic 
wave  fields  to  be  as  shown  in  Figure  1.  The  internal  wave  field  is 
assumed  to  exist  along  a transition  layer  having  a thickness  2L. 


The  internal  wave  field  will  be  of  the  form  £(x,  t)  = b cos  (k  x-u)  t) 

s s 


where  k is  the  horizontal  internal  wave  number  and  ou  /2tt  is  its  angular 
s s 

frequency.  The  perturbation  in  the  index  of  refraction  is  of  the  same 
form,  so  we  shall  write  for  the  temperature  field 


T/6T=  or  cos  (k  x - ou  t) 
s 6 


and  for  the  incident  acoustic  wave 

t (x,z,t)  = i/o  cos  (&.r  - fit), 
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(2) 


(3) 


(4) 


(5) 


vi 


where  K.  r = Kx  sin  9 - Kz  cos  0 and  is  the  amplitude  of  the  incident 
acoustic  velocity  potential  and  f = Q/2tt  is  the  frequency  of  sound. 

We  take  a Fourier  transform  in  time  of  Eqn.  (3)  and  use 

CO 

F^f(x)g(x)|  = j F(y)G(a-y)dy  to  obtain 

- oo 

mO 

A1'tCXJ^£U)+  klit(X,2L,£o)=  (.Xjlj  0>')  Jcj'j 


where  f (iu-uu1)  is  the  Fourier  transform  of  T/6T.  Similarly,  taking 
Fourier  transforms  of  (4)  and  (5)  we  have,  respectively, 


ft"- co')=  (w - 60*  - CJS)  + e.-'ikf(S(cj  + 60'-  60,)  ^ 


dl>(x,T.)(3)=  ^|S'Cco,-X$ecK'r+5'(to'  + IL)e"i:^‘r  j* 

We  will  employ  the  Born  approximation  in  assuming  that  the  scattered 
field  is  not  too  different  from  the  incident  field.  Hence,  we  substitute 
Eqns.  (7)  and  (8)  for  f (tu-tu1)  and  i|r  (x,  z.oo1)  in  the  integrand  of  Eqn.  (6), 
to  obtain 


Eqn.  (9)  shows  that  the  internal  wave  field  modulates  the  incident  acoustic 
radiation  such  that  acoustic  components  with  frequencies  Q ± are 
generated  in  the  scattering  region. 

The  associated  scattered  field  can  be  calculated  using  the  two- 
dimensional  Green's  function 


G (x,ilx.,z.)=  f + (x-zo)*-') 


(10) 


which  is  a solution  to  the  equation 


A'Gr  +kt6  =-SCr-£.) 


(ID 


Eqn.  (10)  represents  the  field  produced  by  a line  source  of  unit  strength 
at  (x0,  z0).  The  scattered  acoustic  field  from  the  internal  wave  is  then 
represented  as 


<o)=  £ 00|z0*  to)  G (x,  z | X0,zo)  d L0  > 

-CO  - L 

where  we  define  the  source  S(x0,  z0  ; u>)  by  the  right  hand  of  Eqn.  (9). 
viz.  , 


(12) 


S C X.,  i U>)  (to-jd-  to,) 
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+ S(<j-SL  + o,)+e±(*-Z-'**-)  + glu+sl.c4>erX(t-?.-kx.)^ 

We  employ  the  Fourier  integral  transform,  Morse  and  Feshbach  (22) 


H.(krU 

4 


i fju  Cj i/  + KzCz-z*)] 

= \<iKAdK2 

J Kl+Kl-k* 


to  replace  the  Hankel  function  in  Eqn.  (12).  The  path  of  integration 
is  chosen  so  as  to  lead  to  an  outgoing  wave  from  the  source  point  r = 0. 
Integrating  Eqn.  (14)  we  obtain 

* //>)=  1 (JK,  S^lL 

and  the  inverse  Fourier  transform  of  this  is  found  to  be 

trJ  • « - • 

The  integral  over  x0  is  performed  first.  From  Eqn.  (13)  we  notice 

that  the  source  term  contains  exponentials  of  the  form  exp  ± i (K  +k  )x 

x s c 

and  therefore  f.  = i (K^  t kj.  Hence,  Eqn.  (16)  for  )r  = K t k yields 
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if  we  let  k sin  0 = K sin  0 + k , then  k2  - (K  sin  0 + k )2  = 

V,. ? s s 

k - k sin  0 = k cos  0 and  (18)  may  be  written  as 

s s 

CO 

fe+(K*  (kr)Jx0=  ^[KXSC^&:+  kcotdUz'*-)] 

J * ° k c«0j 

■00 

In  general  we  can  write  for  the  integral  for  p = ± (K  + k ) 

ao 

L±Z  (Kx  ± ks)Xo  ( !<r)d*o  = *?■--  ei[ikxsi*  9}  + k coiG,*  (z-zj] 

I ° KCoaSj'  ^ 


where 


is  . ± 

~ sin0  = K sin  0 ± k 
c s s 


The  integral  over  z0  is  performed  next.  The  source  term  contributes 

exponentials  of  the  form  exp  (±  iK  z0).  Obtaining  the  companion  exponen* 

z 

tial  from  Eqn.  (19)  the  integrals  under  consideration  are 


L _ 

J [K.-k-tf] 


f*“ 


-i[Kz  + kcos0*Jz.</z,  _ 


IL  ^[Kz + koosB*]L 

[Kj  + kc«i  9+J 


Combining  now  the  results  of  Eqns.  (19),  (21)  and  (22),  Eqn.  (12) 


yields 
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ft  L^i- 


J S (*>-£.-*$ -I ,i£2*sU6!  + 2*‘”<>l]  sUtKt-£c~9t]L 

1 [K,-f-»;]L 

+ S(m-q^>-  ■ JL*'**  +*““*]  j~ [K»-  ^co.gjJL 

““*«■  [K2-^«5e;]L 


+ S^+A+^-i- «' 


+ 8 (cj  +SL-  t^s) — ! e 

7“*®* 


£“*<£  [Kz  + £ cviBtjL 

i -{9^yt^e:+^*cM8;]^rj^+^c«e-]L'| 


[K2  +■“  cm  9 j]  L 


RESULTS 


Taking  now  the  inverse  Fourier  transform  we  have 

DO 

$ (x,  z,  t)  = J j (x,  Z,  0) ) e iUJtduu 

-•e 

which  yields  after  we  perform  the  integration,  using  Eqn.  (23) 


ft 


-Nv 


(A+  cjs)iJ. 


where 


t*- 


LflNi*  w[4-e-^-eyL 

4t*  S=|iiic«e;[fi«1e.S^c«0*jL 
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(23) 


(24) 


(25) 


(26) 


4 


* 

K*-r 


w e*  + t c,  et 


and 


)j . r = z9z±^h.  xsi^ei  - ^ 


Z CoS  &* 


Eqn.  (25)  shows  that  the  scattered  acoustic  field  is  made  up  of  two 
pairs  of  plane  waves  traveling  in  four  different  directions.  The  pair  given 
by  the  first  two  terms  in  Eqn.  (25)  represents  the  forward  scattered  waves 
while  that  given  by  the  second  two  terms  represents  the  backscattered 
waves. 

In  addition,  from  Eqn.  (20)  we  obtain  the  Bragg  formulae.  With 

k = ± (Cl  ± ) we  have 

s 


ou 


± n sin  0 . s 

S1"  es  ' ± n±V 

S S 

for  the  forward  scattered  waves,  and 


U) 


sin  0+  = 
s 


fl  sin  6 
n±  ft) 

s 


n±  id 


s 


for  the  backscattered  waves.  Eqns.  (29)  and  (30)  give  the  direction  of  the 
scattered  waves  while  Eqn.  (25)  gives  their  magnitude. 

From  Eqn.  (26)  one  finds  that  the  amplitude  of  the  scattered  acoustic 
wave,  for  nearly  normal  incidence,  is  proportional  to  the  perturbation 
in  the  acoustic  path  length  (oL)  caused  by  the  temperature  fluctuation. 


(27) 


(28) 


(29) 


(30) 
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DISCUSSION  AND  CONCLUSIONS 


When  an  acoustic  wave  field  crosses  a thin  layer  of  thickness  2L 

perturbed  by  the  presence  of  an  internal  wave  field,  the  acoustic  wave 

will  be  modulated  by  the  internal  wave  and  scattering  and  refraction  will 

occur.  Simple  analysis  shows  that  the  frequency  spectrum  of  the  scattered 

wave  is  resolved  into  a forward  and  a backscattered  component.  Further, 

each  component  is  comprised  of  two  side  frequencies  (Cl  ± u)^)  where  Cl 

and  id  are  the  angular  frequencies  of  sound  and  internal  waves,  re- 
s 

spectively.  For  oceanic  conditions  and  acoustic  radiation  with  carrier 
frequencies  in  the  audible  part  of  the  spectrum,  Doppler  shifts  of  a few 
tens  of  Hertz  are  expected.  Recent  experimental  results  obtained  from 
Project  MIMI  appear  to  indicate  the  existence  of  a Doppler  shift  that  has 
the  same  shift  frequency  as  that  that  would  have  been  produced  by  internal 
waves  in  the  region. 

For  the  case  in  which  both  free  surface  waves  and  internal  waves  are 
present,  the  resulting  spectrum  of  the  forward  (or  back)  scattered  com- 
ponent will  be  comprised  of  the  fundamental  acoustic  frequency  and  four 
additional  side  frequencies;  the  outer  two  resulting  from  surface  gravity 
waves  and  the  inner  two  from  internal  waves. 
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